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CN I Optimal set of 1050 rovibronic energy levels for 35 triplet electronic states of D2 

■ has been obtained by means of a statistical analysis of all available wavenumbers of 
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ly-^ . triplet-triplet rovibronic transitions studied in emission, absorption, laser and anti- 
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crossing spectroscopic experiments of various authors. We used a new method of the 
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analysis (Lavrov, Ryazanov, JETP Letters, 2005), which does not need any a priory 



^ • assumptions concerning the molecular structure being based on only two funda- 

o : 

^ , mental principles: Rydberg-Ritz and maximum likelihood. The method provides 
O 

I the opportunity to obtain the RMS estimates for uncertainties of the experimental 
>>■ 

(-H ' wavenumbers independent from those presented in original papers. 234 from 3822 



published wavenumber values were found to be spurious, while the remaining set of 
the data may be divided into 20 subsets (samples) of uniformly precise data having 
close to normal distributions of random errors within the samples. New experimen- 



> 

■ tal wavenumber values of 125 questionable lines were obtained in the present work 



Optimal values of the rovibronic levels were obtained from the experimental data set 
consisting of 3713 wavenumber values (3588 old and 125 new). The unknown shift 
between levels of ortho- and para- deuterium was found by least squares analysis of 
^ . the a'^S^, V = 0, N = -\r 18 rovibronic levels with odd and even values of N. All 



the energy levels were obtained relative to the lowest vibro-rotational level {v = 0, 
= 0) of the 0^5]+ electronic state, and presented in tabular form together with the 
standard deviations of the empirical determination. New energy level values differ 
significantly from those available in literature. 
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I. INTRODUCTION 

The present work is devoted to the studies of the wavenumbers of the triplet-triplet 
radiative electronic-vibro-rotational (rovibronic) transitions and empirical determination of 
the optimal set of the triplet rovibronic energy levels of the D2 molecule. 

Studies of spectra and structure of molecular deuterium represent not only understand- 
able general interest (an isotopomer of simplest neutral diatomic molecule), but have also 
direct practical value in connection with wide use of D2 in physical experiments and in 
various technical applications: from gas-discharge sources of ultraviolet radiation [l^ up to 
tokamak- reactors 

The spectrum of the D2 molecule is caused by both singlet-singlet and triplet-triplet rovi- 
bronic transitions. The intercombination lines were not observed yet. The most interesting 
resonance singlet band systems are located in vacuum ultraviolet (VUV), while the triplet 
transitions are responsible for major part of light emission of ionized gases and plasma in 



near infrared, visible and near ultravio 



of non-equilibrium plasmas , l6|]. Therefore, the triplet rovibronic levels and 

optical transitions between them were chosen as the object of the present research. 

The energy level values evaluated from experimental wavenumbers are often called ob- 
served energy levels irrespective of a method used for the data processing. Such rather 
disputable terminology implicitly assumes that all spurious experimental results and system- 
atic errors are absent, and random errors are negligibly small (Otherwise, several different 
sets of the observed energy level values may be obtained from the same set of experimental 
wavenumbers by applying different methods of the data processing.). Therefore, we prefer 
to use the term empirical level values ^ (for any values obtained from experimental data) 
and optimal level values (for the values obtained by pure statistical approach [s], 0]), leaving 
the word observed for experimental values of the wavenumbers which are unambiguously 
linked with the wavelengths — real observables in spectroscopy. (It should be also men- 
tioned that in molecular spectroscopy energy levels in cm~^ are often called molecular terms 
(see e.g. jl^). Absolute values of both levels and terms are always positive and calculated 
relative to the lowest vibro-rotational level of the ground electronic state considered as a 



et. They are often used for spectroscopic diagnostics 



^ Compare with the term experimental level values in 
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zero level. That is in contrast to original meanings of terms and energies adopted in atomic 
spectroscopy.) 

Generally speaking, the best set of empirical values of energy may be named a set of 
optimum values (optimal for available experimental data) only if these values do not depend 
on a procedure of their obtaining, and are entirely defined only by amount and quality 
of available experimental data. In such favourable case the energy level values may be 
considered as secondary experimental data. 

It is well known that diatomic hydrogen, being the simplest neutral molecule, has a most 
sophisticated emission spectrum. The hydrogen band spectrum, caused by spontaneous 
emission due to rovibronic transitions, does not show a visible, easily recognizable band 
structure, but has the appearance of a multiline atomic spectra. The peculiarity of molecular 
hydrogen and its isotopic species — abnormally small nuclear masses — leads to high values 
of vibrational and rotational constants and large separation between vibrational and rota- 
tional levels of various excited electronic states. As a result, various rovibronic spectral lines 
belonging to different branches, bands and band systems are located in the same spectral 
regions, leading to the overlap of various band systems, bands and branches, as well as the 
mixing of rovibronic spectral lines having different origins. The small nuclear masses stimu- 
late a breakdown of the Born-Oppenheimer approximation due to electronic- vibrational and 
electronic-rotational perturbations having both regular and irregular character; this combi- 
nation seriously complicates the interpretation of the spectra of hydrogen isotopomers and 
the unambiguous identification of rovibronic spectral lines. Symmetry rules for permutation 
of identical nuclei in homonuclear isotopomers (if 2, -D2 and T2) cause the known effect of 
the intensity alternation of neighbouring lines within the rotational structure of bands due 
to the alternation in degeneracy of successive rotational levels with odd and even values of 
rotational quantum number (e.g. 1:2 in the case of 1^2) • This effect also masks the visible 
structure of branches resulting in serious additional difficulties for identification of rovibronic 
spectral lines. 

Thus, most of the lines in the optical spectra of hydrogen isotopomers have not yet 
3een assigned in spite of tremendous efforts by spectroscopists over the previous century 

As an example, in the latest compilation of experimental data 
for molecular deuterium D2 [3] , the working list of 27488 recorded lines contains only 8243 
assignments. These assignments were obtained by traditional methods of analysis using 
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wavenumber combination differences (method of common differences) and Dunham series 



expansions 
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13| . sometimes together with comparison of molecular constants obtained 



for different isotopic species p^. Later on, the traditional methods were supplemented by 
taking into account the line intensities 18| and the results of ab initio and semi-empiric 



calculations [161, 1171 
Recently in [19| 



the identification of rovibronic lines of three singlet VUV band sys- 



tems of D2 has been obtained by common use of semi-empiric calcu 
the transition probabilities) together with the program IDEN 21 



ations (of the levels and 



22| developed for visual 



identification of complex atomic spectra. Although the approach seems to be a rather power- 
ful and prospective tool for finding new assignments in the spectra of hydrogen isotopomers, 



the final results of [19|, l20|] can not be considered as optimal values o: 



the rovibronic energy 



levels for the -D^II^, Z^'^II^,, and B'^J2u electronic states studied in [l9|, |20]. They were simply 

calculated as the average of values derived from different observed line wavenumbers added 

;o the corresponding lower vibro-rotational levels of the X^T,^ ground electronic state from 

14| . assumed to be known with negligible uncertainty. Actually, those lower levels were 

obtained from wavelengths of two other VUV band systems {B^T,^ —X-^H^, C^H^ —X-^Ti^) 

by two independent methods of the analysis which gave noticeably different results (compare 

columns 2 and 4 of table 4 from [l^). It is evident that common use of experimental data 

for five singlet band systems mentioned above may give another set of rovibronic levels for 

X^Tif ground electronic state. Then the level values for the exited states will be changed 

— I I — I 

somehow. In our opinion, the assignments found in 19|, l20| also should be confirmed by 
the appropriate statistical analysis of all available in literature experimental data on the 
wavenumbers for all other singlet band systems of Do. For that purpose it is possible to use 



23 



24|. 



the approach proposed and realized in our papers 

Currently, available information concerning triplet rovibronic energy level values of D2 
molecule exists in the form of the list of molecular constants for Dunham series expansions 
in 
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25[ |. and tables of rovibronic levels obtained in [l6l |. 



In the case of the hydrogen molecule the Dunham coefficients are known to provide 
a very poor description of the rovibronic energy level values 26|]. Moreover any set of 



molecular constants depends not only on molecular properties but also on the theoretical 
model used in the analysis of the spectrum (the type and length of the series expansion, form 
of a model Hamiltonian, etc.), as well as on amount of measured wavenumbers and even 
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on the distribution of the spectral hnes over the vibrational and rotational levels (see e.g. 



BE 



271]). Therefore, the molecular constants are rather conditional and the procedures 



used for their determination are subjective in character. In practical use of reference data on 
molecular constants, the question arises as to how many rovibronic levels of a particular type 
can be described by a given set of constants with a required accuracy. In the best case, this 
question can be answered only by the researcher who obtained the given set of molecular 
constants. However, series expansions and molecular constants may be useful both in a 
process of classifying rovibronic lines in unknown spectra and in calculations which do not 
pretend for a high accuracy. 

The data reported in 16|] , in general, give a rather good description of the D9, spectrum, 
but they are also not free from criticism. The method of the analysis used in IGj is based 
on the common use of the combination differences, some selected wavenumbers for certain 
transitions and one by one multistage treating of separate branches, bands and band systems. 
The sequence of the steps chosen in [l6| is not the only possible analytical arrangement. 
Therefore, the data thus obtained can not be considered as an optimal set of levels providing 
the best description of all observed wavenumber values. On our opinion the reported in |l6l | 
uncertainty a = 0.05 cm~^ for all levels and all wavenumbers of lines is the overestimation 
caused by the randomization of various deviations when working with many thousands of 
levels and lines. It should also be mentioned that after publication of 16| , new experimental 
data on the wavenumbers appeared 28l. l29|. 

Recently, a new method of a pure statistical analysis of experimental data on the rovi- 
bronic transition wavenumbers for empirical determination of the optimal energy-level values 
has been proposed [8|] and successfully applied for the derivation of rovibronic level values of 



the singlet states of the BH 



231] and the triplet states of the H2 2J] molecules. The method 



is based on only two fundamental principles: Rydberg-Ritz and maximum likelihood. This 
approach differs from known techniques in several aspects: 1) does not need any assumptions 
concerning an internal structure of a molecule; 2) does not employ intermediate parameters, 
such as molecular constants in the traditional methods; 3) a one-stage optimization proce- 
dure can be used for all available experimental data obtained for various band systems, by 
various methods and authors, and in various works (Simultaneous adjustment of all levels 
eliminates the possibility of accumulative errors in traditional methods, caused by multi-step 
treating of various bands one by one, as well as by transfer of any occasional error for one 
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level, appeared on one step, for all higher levels obtained during further steps); 4) provides 
the opportunity of a rational censoring of the experimental data in an interactive mode (thus 
allowing the user the option to eliminate obvious errors (outliers) and misprints, to revise 
incorrect line assignments, and to compare various sets of experimental data for mutual con- 
sistency); 5) gives an opportunity of independent estimation of experimental uncertainties 
by analyzing the shapes of error distributions within various samples of the experimental 
data (It means that our statistical estimates of the root-mean-square (RMS) experimental 
errors are totally independent on the estimates reported in original experimental papers); 
6) provides an optimal set of rovibronic level values as well as the uncertainties of their de- 
termination (standard deviations SD and the covariance matrix) caused only by the amount 



and quality of existing experimental data 

A necessary precondition for an application of the method [8| is an existence of more or less 
dependable assignments for majority of rovibronic spectral lines under the study. Therefore, 
the new method should be considered as complimentary for the traditional and new methods 
of identifying unknown spectra. However, since the new, pure statistical method is objective, 
it is better to consider any identification of a rovibronic spectrum as preliminary until it is 
not corroborated by that, the use of certain assignments in the framework of this method 
does not lead to contradictions with the Rydberg-Ritz principle or with normal distribution 
of random experimental errors. 

The goal of the present paper is to report the results of applying the new method 81] for 
statistical analysis of the rovibronic spectral line wavenumbers of triplet band systems and 
determining the optimal set of rovibronic energy levels for all known 35 triplet electronic 
states of molecular deuterium: a^S^ , c^n+, c^H", rf^n+, (i^II^, e^S+, /^S^l;, 5'^^+, /i^S^J;, 
t^U+, i^U-, fA+, j3A;, Pn+, fc^n-, n^Ut, n^n-, q3j:+, r^U+, r^U-, s'A+, s^A;, 

u^Ut, u^U-, {rpfUt, {rpfU-, {SpfUt, {SpyU-, {QpfUt, {dpfU-, (6rf)3S+, (7rf)3S+, 
(8d)3E+ and (M)3E+. 

II. STATISTICAL APPROACH TO EMPIRICAL DETERMINATION OF 
OPTIMAL VALUES OF THE ROVIBRONIC ENERGY LEVELS 



In principle an analysis of experimental data on wavenumbers of spectral lines may be 
considered as consisting from two separate parts: 1) an identification (classification, assign- 
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ment) of spectral lines (an establishment of a correspondence between observed spectral 
lines and pairs of the initial and final levels responsible for occurrence of these lines in a 
spectrum), and 2) the determination of energy level values from measured wavenumbers 
(or some other magnitudes related to the levels, like molecular constants 26|], or potential 
curves SOj). The aim of the first stage of an analysis is to find most likely identification and 
to prove correctness of the assignments. The goal of the second is to obtain the best possi- 
ble values of the energy levels as well as the uncertainties of their empirical determination. 
In spite of the difference in objectives these two parts of an analysis are closely connected 
because quality of final results in both cases is determined by their ability to reproduce 
observed spectra. Therefore, quite often the parts are not distinguished in original papers. 
The main part of the present work represents the second stage of the analysis of the 



;ri 



161, 



et-triplet spectrum of Dp, molecule, based on the results of the first stage reported in 



28 



29 



31 



32 



34 



35|. ^ Therefore in the further consideration it is supposed, that 



more or less reliable identification is already established for the majority of lines. Then not 
numerous cases of wrong assignments may be revealed as outliers and corrected during the 
statistical analysis (see below). 

All spectroscopic methods of empirical determination of energy levels of atoms and 
molecules are based on the Rydberg-Ritz combination principle corresponding to the Bohr 
frequency rule in quantum mechanics. For rovibronic transitions of diatomic molecules be- 
tween electronic states corresponding to Hund's case "b" for angular momenta coupling 
(negligible multiplet splitting) it may be written as: 



^n'V^N" — En'v'N' — En"v"N", (l) 

where n indicates an electronic state, and v — the vibrational quantum number. In case of 
diatomic hydrogen isotopomers the multiplet splitting and hyperfine structure of levels and 
lines are negligibly small and usually unresolved in experiments, therefore in the present 
work rotational levels are characterized by the quantum number of the total angular 
momentum of a molecule excluding spins of electrons and nuclei. Upper and lower rovibronic 
levels are marked by single and double primes, respectively. The u^i'F^l!^^,, is the wavenumber 



^ In some cases discussed below, we have had to do the job of the first stage when it was necessary to change 
earher assignments and to obtain new experimental values of the wavenumbers of some questionable lines. 
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(recalculated for vacuum conditions) of the n', v', N' - n", v", N" rovibronic transition, and 
EnvN are corresponding energy level values. 

One may see that the experimentally measured wavenumbers of the spectral lines are 
related only to the differences between pairs of rovibronic energy levels. So the values of 
the levels cannot be directly obtained by means of the Rydberg-Ritz principle only. There- 
fore, traditional methods of empirical determining the rovibronic energy levels require the 
introduction of some additional (with respect to the Rydberg-Ritz combination principle) 
assumptions regarding properties of molecules and, hence, are semiempirical and allow some 
subjectivity of a researcher (see details in j^). 

The classification of the levels and lines by certain sets of quantum numbers (n, v, N) 
is important for the assignment of the wavenumbers to the rovibronic levels between which 
the transitions occur and, in particular, for the application of the selection rules. When this 
assignment is carried out, the specific designations of the levels and lines are not important 
any more. The rovibronic levels may be denoted by natural numbers in an arbitrary order. 
The notation EnvN may be changed to a more compact notation Ei, and the wavenumbers 
are designated by pairs of indices corresponding to the initial and final levels. Then the 
expression ([1]) becomes: 



E,-E,. (2) 



Suppose that the set of available experimental data consists of n,, wavenumber values 
obtained for Ut transitions caused by the combinations of Ue energy levels, riy > tit > 
because an amount of the transitions allowed by selection rules usually considerably exceeds 
quantity of combining levels. Furthermore, experimental data for the same transitions may 
be obtained and reported in various publications. One may insert all experimental wavenum- 
bers into left hand side of the ([T]) one by one, and consider the energy levels as adjustable 
parameters to be obtained from the experimental data. Then ([T]) turns into a system of Uy 
equations with ue unknown quantities E^. This system of equations is overdetermined and, 
hence, is inconsistent because the experimental data always involve measurement errors. 

Straightforward general solution for solving such problems is well-known in mathematical 
statistics. That is the least squares method based on the assumption of finite second moment 
of the distribution function for random errors 3^. In our case it consists in the minimization 
of the weighted mean-square deviation between observed v^J^^ and calculated (as differences 
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of adjustable energy levels Ei, Ej) values of rovibronic line wavenumbers, or the sum 



E 



{E,-E,)- 



(J expt 



(3) 



The values a ^^vt are the RMS estimates of experimental errors (one standard deviation 
— SD) for each experimental datum, and the summation is performed over all available 
experimental data. Due to the linearity of the equations ([2]), the optimization problem 
comes to solving a system of linear algebraic equations. The complete solution involves the 
inversion of the nE*nE sparse matrices. When the number of desired levels is of the order 
of a thousand, as is usually the case, this problem can be solved even using modern personal 
computers, (see [9] for details)). If the experimental errors are random and distributed 
according to a normal (Gaussian) law, the obtained solution corresponds to the maximum 
likelihood principle 36|. 

The attempts to apply such statistical approach for determination of the energy level 
values appeared almost simultaneously in atomic (e. g. the spectra of the Er I, Er II 



37l |. CI II [38| by the inversion method, and the spectra of the Si I |39|], CI I |40|] by the 
iterative method) and in molecular |4l| spectroscopy when first digital computers became 
available. It is interesting to note, that in spite of obvious disadvantages of the iterative 
method (the problems of convergence and lack of the covariance matrix (see p. 6 in 3]), its 
usage can be met in current publications (see e.g. |43|). 

In several studies reviewed by Aslund in 41|, it was proposed to determine the rovibronic 
levels of molecules on the basis of the combination principle by solving the overdetermined 
system of equations with the use of a computer. Owing to the limited capabilities of the 
computers of the day, the optimization procedure was realized for processing the wavenum- 
bers for separate bands only. As a result, those authors overlooked the opportunity to 
overcome the problem of the limitations for the method of common differences caused by 
the Laporte's selection rule allowing rovibronic transitions only between levels with different 
+ and - parities [sl. (Only recently, it was shown that, if it is possible to use experimental 
data on rovibronic lines that pairwise couple three or more different electronic-vibrational 
states, then the system of equations ([T]) contains the link between the line wavenumbers and 
the values of all rovibronic levels involved [8]. Thus, in the case of heteronuclear molecules 
the above problem of the existence of uncoupled sets of the levels with even and odd rota- 
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tional quantum numbers disappears.) Consequently, the relative positions of unrelated odd 
and even rotational levels, as in traditional approaches, were obtained using the Dunham 
approximation and molecular constants. It should be noted that the attempts of direct 
empirical determining the rovibronic levels in 4l[ were inconsistent from the very outset. 
The level values were considered as some intermediate parameters of a molecule, and the 
procedure of their determination was only one of the steps on the long way from the mea- 
surement of the line positions on a photographic plate up to obtaining a particular set of 
vibrational- rotational constants for various electronic states. Later on, Aslund 44 1 actually 



abandoned his initial idea in favour of so-called direct approach [27] based on the one-step 
optimization of sets of vibro-rotational constants for upper and lower electronic states of a 
bad system. 

The empirical determination of the energy level values by means of the least squares 
method looks natural and simple, but its correct realization in practice encounters serious 
difficulties. The matter is that the result of the minimization of (IHl) (i. e. the optimal 
set of the level values and the matrix of co-variances) depends not only on the values of 
experimental data, but also on weighting of various data (i. e. from what values of the error 
estimates cr^exvt are included into the input data set). Therefore, for the correct optimization 
of the adjusted parameters (required level values) it is necessary to know a dependable 
RMS error estimate for each experimental wavenumber value, or to have the set of the 
wavenumbers consisting from several subsets (statistical samples) of uniformly precise data 
with known RMS error estimates for each sample. 

Unfortunately, authors of experimental works usually limit themselves to some remarks 
of general character (concerning resolving power of a spectrograph, linear dispertion, typical 
line widths, etc.) and, at the best, to some rough estimates of an order of magnitude or 
the upper limit of possible experimental uncertainty. Most often the real accuracy of the 
reported data is uncertain. The situation in the literature containing experimental data 
analyzed in the present work will be illustrated below, but it should be mentioned that the 
same uncertainty may be met even in current publications. (To not be unsubstantiated, very 



recent paper 



431] may be used as a typical example. In the part describing experimental 



setup one may read: "The estimated error is -1-0.005 A for single lines, but many lines are 
blended in complex emission peaks."; while in the part describing data processing one may 
find: "As input to the code, 1314 classified lines were used with uncertainties on their wave 
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numbers smoothly decreasing from 0.33 to 0.10 cm" between 1150 and 2800 A." Attentive 
reader can easily recognize that the wavenumber uncertainties correspond to wavelength 
uncertainties 0.0044 and 0.0078 A not equal to those for single lines. But a lot of questions 
remain. What about blended lines? Why the weighting of various data depends only on 
wavenumber values? What kind of the smooth decrease was chosen and why? We hope 
that the authors have certain answers and did their best to get best possible level values 
they are interested in. But what an independent researcher can be able to do with those 
wavenumbers after several decades?) 

The problem of weighting various experimental data according to their accuracy in the 
framework of the least squares fitting was mentioned from the very beginning [37]. But in 
practical applications it was always considered as an issue of secondary importance, leaving 
a lot of room for author's subjectivity. May be it is because researchers use to process 
their own experimental data and have certain opinions concerning their accuracy. Then, 
irrespective of the way of the data processing, the empirical level values are considered 
as optimal when they are able to reproduce experimental wavenumbers with an accuracy 
corresponding to author's estimates for the experimental errors ("chosen tolerance"). That 
is certainly reasonable, but only if author's expectations coincide with real experimental 
uncertainty caused by all possible sources of errors. Thus, in spite of principle objectivity 
of the least squares method the results of its application have subjective character reflecting 
author's individuality. Achilles' heel of all known to authors of the present work papers 
devoted to the determination of the level values by means of the least squares adjustment 
is ignoring or subjective estimating of the experimental errors. 

Development of pure statistical approach to the problem of empirical determination of 
energy level values was motivated by author's desire to overcome the subjectivity dis- 
cussed above. Main ideas of the approach may be formulated in the following way. Since 
optimal values of rovibronic levels are to be entirely defined only by amount and qual- 
ity of available experimental data, a procedure of their obtaining should be based only on 
the fundamental principles. Rydberg-Ritz combination principle provides the link between 
experimental wavenumbers of rovibronic transitions and desired energy levels. Maximum 
likelihood principle gives the optimization criterion in the form ([3]) for the normal distribu- 
tion function of random errors. The set of input data for minimizing ([3]) should be prepared 
by means of statistical analysis: 1) The data having systematic errors and spurious results 



12 



(outliers) should be uncovered and eliminated; 2) The whole set of available experimental 
data should be divided into finite number of subsets (samples) of uniformly precise data 
with normal error distributions; 3) the values of an experimental RMS uncertainty for the 
samples (unique for all data belonging to the same sample) are obtained, being independent 
from those named in original papers. Optimal set of the rovibronic level values is obtained 
from the prepared input data set in the one-stage optimization procedure. The accuracy of 
empirical determination of the optimal level values is characterized by the covariance matrix 
D{E), square roots of its diagonal elements are used as the uncertainties (one SD) of the 



level values {cte, = \J Du^E)). (The method may be easily generalized for other types of the 
distribution function of random errors.) 

The computer code developed for practical realization of the method is based on the 
minimization of the weighted mean-square deviation ([3]) by the inversion of the nE*nE sparse 
matrices . The z/*^f^* and a^^expt values are the input data, while the optimal set of the 
energy levels Ei and their covariance matrix D{E) are the output data of the minimization. 
The later is used for calculating and visualization of various quantities suitable for interactive 
statistical analysis of the experimental data including rational censoring of various groups 
(samples) of the data and the derivation of the RMS uncertainties of experimental data 
independent from those reported in original papers. 

The interactive analysis consists in the studies of the shape of the distribution functions 
of the weighted unbiased deviations 



for various sets (statistical samples) of the experimental wavenumbers, k being the number 
of the sample. Here the u^^^'^ is certain experimental datum for the wavenumber of the i 
- j transition, while the level values E^ and Ej are obtained by the minimization of ([3]) 
without using the particular datum under the consideration. The cXj^fc is adjustable value of 
the estimate for experimental uncertainty (RMS error estimate) common for all the data 
included into the sample. The dimensionless deviation characterizes a degree of co- 
ordination of the particular experimental datum u^^^^ with all other experimental data in 
the framework of the Rydberg-Ritz principle. Among all, the most important are certainly 
the data obtained for the same rovibronic transition in various works, and the wavenumbers 
of other lines directly connected with the upper n', v', N' or the lower n", v", N" rovibronic 






(4) 
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levels. 

The values Ei and Ej may be expressed through the corresponding level values and the 
covariance matrix D[E) obtained with the use of all experimental data, then 



^'''tt-iD,,{E) + D„{E)-2D,j(E)) 
where Dij{E) are elements of the covariance matrix. This makes it possible to obviate the 
need for the repeated minimizing ([3]). 

The empirical accumulative distribution function for each data sample is defined by 

m = ^j:mj<o, (6) 

where rii, is the number of elements in the sample, I{A) — the indicator of event A.^ Here 
the sum is over all data included into the sample. After each minimization of (3) it is 
possible to calculate and to see F{^) for any sample of the ^ij values chosen by the user 
of the computer code. That gives an opportunity to study various data samples selected 
from the complete set of all experimental data by some physical considerations: the data 
reported in the certain paper, wavenumbers of lines from the certain interval of wavelengths, 
or belonging to the certain band, or to the certain band system. 

It is evident that if all the data i'^^^'^ included into some sample are free from systematic 
errors and have normal distribution of random errors with RMS estimate equal to cr^k, then 
the empirical distribution F{C,) should be close to normal cumulative distribution function 
with the zero mean and variance equal to unity, namely 



1 2 

m) = ^ / expi-^)dx. (7) 



— oo 

Therefore, it is useful and convenient to provide the interactive analysis of experimental 
wavenumber values by studies of the shape of the empirical functions 

m = Fo\F{0), (8) 



' Ii^^J < = 1 for < ? and 1^,, < - for > ^ 
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calculated for various groups (samples) of available data. Here Fq ^ is the function reverse 
toFo(0,i.e. F,-\Fo{0)^C 

If the empirical distribution function F{C,) is close to normal distribution, the empirical 
function (IH]) is close to linear function 



m = e 



(9) 



The way of dividing (sampling) the whole set of experimental data onto the subsets is 
not determined leaving some subjectivity for a researcher. But the shape of the function 
$(^) is objective characteristic of the sample being determined by experimental values of 
the wavenumbers, Rydberg-Ritz principal and only one arbitrary parameter 0",^^,. If for a 
certain data sample adjusting the value of the a^k leads to the shape of <l>(^) fulfilling the 
criterion ([9]), the data included into the sample are to be considered as measured with 
uniform precision and characterized by the normal distribution of random experimental 
errors. The obtained value aj^k is their experimental uncertainty (RMS error estimate). 
Thus, the method of a posteriory statistical analysis proposed in 8|, |9| gives the objective 
mean to get proper results by the studies of the $(0 functions, trying various kinds of 
sampling and adopting or rejecting various guesses. First applications of the approach in 
our recent studies of rovibronic spectra of BH and H2 molecules 8|, l9|, l23|, showed that the 
method allows revealing and rejecting both single outliers (caused by wrong assignments, 
blending, spurious readings and misprints) and large data sets, including data__of some 
experiments (probably caused by systematic error of the wavelengths calibration 



III. STATISTICAL ANALYSIS OF PUBLISHED WAVENUMBERS OF 
TRIPLET ROVIBRONIC TRANSITIONS 



All available values of rovibronic transition wavenumbers studied in emission. 



laser and anticrossing spectroscopic experiments of various authors 
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28 



29 



absorp 



31 



;ion, 



32 



34 



33, 



35| were analyzed in the present work. In overwhelming majority of the works the fine 



and hyperfine structures of lines have not been resolved, and the reported experimental 
wavenumbers correspond to intensity maxima of the observable line profiles. In cases of 



partly resolved structure of lines |28|, l29| we used wavenumbers of the brightest components. 



The all identified triplet electronic transitions are naturally divided into three groups of 
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the band systems. These are: i) the 9 band systems having one common low electronic 
state {n^Ag — ct'^S^, with A = 0, 1 and n = 3 — 9); ii) the 12 band systems having another 
common low state (n^A^ — c^n„, with A = — 2 and n = 3 — 9); and iii) the 7 band 
systems connecting various n^Ag and n^A„ electronic states (with A = — 2 and n = 2 — 4) 
including the — c^Hu transitions. Therefore, the whole set of available experimental 

data possesses different informational content about rovibronic levels of different electronic 
states. Thus, the information on rovibronic levels of the c^Hu electronic state is contained in 
the wavenumber values of the 13 band systems, about the o^S^ state — 10 band systems, 
and about for example the e^T,^ state — 6 band systems. For some electronic states as a 
source of the information about rovibronic levels we can use wavenumbers of lines from very 
few or even only one band system. In such cases the obtained values of the rovibronic levels 
are essentially less reliable. 

Most of the data used in our statistical analysis are those collected from earlier works 
(partly unpublished) and reported in |l^. This compilation contains wavenumbers of 3117 
rovibronic spectral lines assigned as triplet-to-triplet transitions. 31 from them have unas- 
signed upper electronic and vibrational states^ They were not used in the analysis, so 
only 3086 experimental data were taken from [IGj]. The 83 wavenumber values obtained in 



33| coincide with those reported in 



16| (although 12 of them appeared in Appendix C of 



161] as unassigned). To prevent doubling of the same experimental results we excluded 71 
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32| and [l6| are obtained 



wavenumbers from the data sets of [33] . The data reported in 
for the lines of the same band systems, but wavenumber values are different. Taking into 
account random character of the differences we decided to consider the data of those papers 
as the results of independent experiments. Thus, the initial data set used for the start of the 
statistical analysis contains: 350 wavenumbers from paper 



from 1341, from 



, 285 from jSlLall 37 data 
33|, 1 from [35], 3086 horn Q, 81 from [281, and ^froniQ . It should 



be stressed that the wavenumbers of spectral lines used in the present work are spread over 
the very wide range 0.896 — 28166.84 cm~^ from radio frequencies up to the ultraviolet. 

On the first iteration of our statistical analysis we assumed, that all sets of experimental 
data presented in various papers represent the samples of uniformly precise data with RMS 
error uncertainty equal to the estimates presented in original experimental works. In the 
review |l^ there are two remarks concerning accuracy of the wavenumbers. In experimental 
part it is written: "... the widths of the lines themselves limited the attainable accuracy. 
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which was a few hundreds of a cm~^". At the end of the analysis the authors of [l^ 
mentioned: "The mean (0-C) value for all triplets is —0.001 ±0.051 cm~^", and came to the 
conclusion that " The precision of the measurements appear to be roughly uniform from the 
infrared to the ultraviolet". In spite of doubtfulness of such estimations we have accepted 
the value = 0.05 cm~^ in the first iteration of our analysis for all wavenumbers taken 



from [l6[]. In papers 3l|, l32|, experimental errors are not mentioned at all. Therefore, we 



used the same value 0.05 cm ^ taking into account the estimates made by the same scientific 



group. In 33j it is written: "... the errors in the wave number of good lines should not exceed 
a few times 0.01 cm~^". We certainly don't know which lines are good or bad, therefore 



in this case we again took ai, = 0.05 cm~^. For the data taken from 28|, l29|, l35|] we used 
author's estimates = 0.01,0.003,0.0003 cm~\ respectively. Under these assumptions we 
obtained a set of energy levels by minimization of ([3]) and the empirical function $(.^) shown 
in figure [H 

One may see that the distribution function F{C,) is far from the normal distribution -Fo(^) 
($(^) is too far from the dotted straight line representing the case of normal distribution). 
Too many experimental wavenumber values significantly exceed those calculated as differ- 
ences of corresponding energy level values. The distribution is not close to normal even in 
the area of small deviations < 3). We suppose that this disagreement is caused by the 
deviations of real experimental wavenumber errors and their estimates, reported in original 
works. (On this stage of our analysis it was observed that 33 lines with wavenumbers mea- 



sured in jl6] are the combinations of 55 rovibronic levels representing the blocks of levels 
totally disconnected with all other triplet rovibronic levels of D2. Therefore this wavenum- 
bers were excluded from further consideration as useless. Therefore, only 3053 from the 



data reported in 16| are useful for further analysis. 

Thus there is the necessity of independent estimating of experimental uncertainties of the 
measured wavenumbers. To carry out this estimation we need to divide the data set used in 
the first iteration into some groups (subsets, samples) so, that all experimental data inside 
one subset can be considered as uniformly precise. We tried to sort experimental data to 
subsets taking into account various criteria, selecting the data obtained in various papers, 
in various wavelength regions and in various band systems. But acceptable results for error 
estimates were obtained only when we considered each band as a group of uniformly precise 
data. May be that way of sampling is fruitful because the lines of each band are located 
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in narrow parts of the spectrum, and different bands are usually measured and treated by 
different people even within the same scientific group. It should be also mentioned, that 



such approach gave reasonable results in similar analysis of BH [23'] and H2 2J] spectra. 

Each m-th iteration consists of: 1) the minimization of ([3]) with certain m-th set of the 
input data prepared at the last stage of the previous (m-l)-th iteration. It includes certain 
amount of experimental wavenumbers divided into finite number of samples characterized 
by certain RMS estimates for experimental errors in each k-th sample crl^'^^', 2) the output 
data (optimal energy levels E^"^^ and the covariance matrix D^J^^) are used for calculating 
the $''™'''(^) functions for each sample. That gives us an opportunity to analyze the shapes 
of the fI"^\^) distribution function for all the samples and to find outliers. The data with 
^1™'* larger than 3 were qualified as outliers. In all cases this criterion of censoring was less 



strong then commonly adopted Chauvenet's and Peirce's criterions [45j; 3) The new set 
of input data for the next (m+l)-th iteration is prepared by excluding the outliers from 
current (m)-th data set, and by obtaining for each sample new values of RMS estimates for 
experimental errors cxf,™'' by multiplying the af,™ by the factor which should move the a^k 
closer to unity on the next iteration. Usually we used the expression cr^™'' = ^^'^'(k\ 
which is not strict due to the nonlinear link between erf,™-* and cr^™^^\ 

The square roots of the second moments o"^™'' for empirical distribution functions Flf^\^) 
(for each sample in each iteration) are calculated as 



T.i^^-<^^>yFir^idf), (10) 

where < ^Ij^^ >= J2 Fl:^\c,ij^^), the summing is performed over all the data in- 

eluded into the k-th sample. 

The method of independent error estimation of wavenumbers belonging to the same k-th 
sample may be illustrated by considering one typical example — analysis of the (0 — 1) band 
of the e^S+ — a^T^t electronic transition. 30 spectral line wavenumber values for R- and 

a y 



P-branch lines from this band are reported in [16|]. The lines -R18 and R19 are caused by 
the transitions between the (i?e,o,i95 -^a.i.is) and (i?e,o,2o, -£'0,1,19) pairs of rovibronic levels, 
respectively. Mutual positions of these levels can be derived using wavenumber values of 
these lines only, because there are no other lines, which provide information about these 
levels. Then the value (j4]) for these lines is equal to zero. We don't take into account such 
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lines when calculating empirical distribution functions $[,"^''(^). 

The empirical function obtained during first iteration (based on the sample of 

remaining 28 experimental wavenumbers characterized by common estimate o"^;? — 0.05 
cm^^ mentioned in 16|) is presented in figure EJ^a). One may see, that the distribution 
fI^\C,) is far from the normal distribution, which is represented by a dotted line on the figure. 
The square root of the second moment of this distribution a^^f} = 0.56 ± 0.14 = 0.56(14) is 

(2) 

not close to unity. We may try to shift the a^jj to unity by changing error estimation in the 
following way: a^J = af,°V^^^ = 0.028(8) cm'^. 

In the second iteration for the band under the consideration we used the same set of 
experimental data but new value of a^^^ instead of crf,°\ The empirical function ^^^^ is 
shown in figure [2](b). One may see that now the majority of lines shows the deviations 
sufficiently closer to the normal distribution, and the square root of the second moment 
cr^l'' = 0.9(3) is closer to unity. But one point located in the bottom left corner of the 
figure is in contradiction with rest of the distribution. This point represents the ^ value 
for the line PIO. The measured wavenumber value of this line fa'i,'io = 9292.16 cm~^ 
deviates from the difference of corresponding energy level values derived in this iteration 
^a'i,'io ~ (-^eA9 ~ -^M.io) = —0.0950 cm~^. The deviation exceeds our estimate for RMS 
error cr^^'' = 0.028(8) cm^^ more than three times. Thus the experimental wavenumber for 
the PIO line should be excluded from the input data set as the outlier. 

(2) 

In the third iteration we used 27 wavenumber values and common uncertainty a^f! = 
'^uk = 0.028(8) cm~^ The function $i^''(0 thus obtained is shown in figure [2](c). One 
may see that this distribution function is not close to normal distribution function mainly 

.(3) _ ^(2)^(3) 



because = 0.60(8) is too far from unity. Therefore we have to correct the experimental 



uncertainty as a^;: = o"^,^ cr^^\ After certain number of iterations we are coming to a^k = 
0.012(1) cm"^. Then on the next iteration we have empirical ^k{0 function shown in 
figure [2]^d) and the square root of the second moment a^k = 1.0(1). One may see that now 
the experimental distribution function -Ffc(0 is close to the normal distribution function with 
zero mean and the variance equal to unity. 

Thus, the sample of the 27 wavenumbers under the study may be considered as a group 
of uniformly precise experimental data with RMS estimate for random errors a^k = 0.012(1) 
cm~^ It should be underlined that our estimate is less than a quarter of the value 0.05 
cm~^ declared in 1^. 
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In that way we estimated RMS uncertainties of experimental data by the rational cen- 
soring of experimental data and determination of proper values of cXuk- 

After several dozens of iterations described above we came to the conclusion that 234 
experimental data representing 228 spectral lines have to be exclude as outliers. Remain- 
ing 3588 wavenumbers were organized as 317 samples of uniformly precise data with RMS 
uncertainties obtained by the statistical analysis. RMS errors for many bands are so close 
that it seams reasonable to provide the enlargement of samples by merging studied samples 
having close values of RMS error. This also helps to improve statistics because some bands 
have small amount of data. It is reasonable to unite into the same sample wavenumbers of 
various bands with RMS errors getting to the certain wavenumber interval. The amount 
and widths of such intervals are uncertain. The choice of these parameters may lead to some 
subjectivity. On one hand, the widths of intervals should be large enough for the amount of 
data included in each interval be sufficient for statistics. But on the other hand, it should 
be as small as possible, not to disturb the close to normal error distributions obtained for 
separate bands. Our computational experiments show, that the intervals (0.00^0.01) cm~^, 
(0.01 -7- 0.02) cm^^, (0.02 0.03) cm^^, etc. fullfil the requirements mentioned above. It is 
important, that all energy level values, obtained by the minimization of ([3]) with experimen- 
tal data divided into the such enlarged samples of this size, have been changed significantly 
less then one SD of their determination with the original 317 samples. 

After that, the whole set of the 3588 experimental data was reorganized into the 17 
enlarged samples of uniformly precise data from 16|, |28|, l3l|, l32|, |33|, |3^. The RMS error 
estimates for each of these 17 subsets were derived by the method similar to described above. 
The results of the distribution of the data among the 17 samples together with the list of 
RMS error estimates for the wavenumbers of lines belonging to various vibronic bands is 
presented in table 1 of [46]. Also 1 wavenumber from 35] and 3 wavenumbers from 29] 
were used in further analysis as two separate samples (18-th and 19-th) with RMS errors 
reported in original papers, because the errors were unchanged during the analysis. 

The amount of experimental data Ui, in each of the new samples is presented in figure [3l 
It is seen, that the amount of data in each sample obtained by the enlargement is sufficient 
for carrying out the statistical analysis. Moreover, the distribution of data over the error 
steps looks quite plausible: the number of too precise data {a^k = (0.00 ^ 0.01) cm~^) is 
small (~ 9%), most of the data (~ 85%) are measured with realistic accuracy (0.01 -r- 0.06) 
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cm~^, and the quantity of data monotonously falls down with further increase in the error 
estimate a^k- 

Empirical functions $(0 fo^' fi^^st six united samples representing vast majority of useable 
experimental data are shown in figure HI Correspondent values of RMS error estimates for 
each sample are shown in the figure. One may see that the error distribution functions are 
in much better accordance with normal distribution function then those obtained for the 
data samples representing separate bands (see figure [2]). 

Thus, splitting of experimental data into the samples and the subsequent enlargement 
of the samples as a whole were justified. Taking into account high enough statistics of 
the samples we are coming to the conclusion that obtained error estimates with a high 
probability are close to the values of the real experimental uncertainty. 



IV. EXPERIMENTAL DETERMINATION OF WAVENUMBERS FOR 

QUESTIONABLE LINES 

The statistical analysis of the 3822 known to authors published wavenumber values has 
shown that 3588 of them are useable for empirical determination of the level values, while the 
234 experimental data (concerning 228 spectral lines) have to be excluded as outliers. 161 



of those questionable lines are located within the wave 
for us. 26 from them are the data reported in old 



ength range 4300 ^7300 A available 



31 



321 1 papers and the wavenumber 

n 

values for these lines were measured once more and corrected in the compilation [16|. For 
remaining 135 lines we decided to provide independent experimental determination of their 
wavenumbers. For that purpose we use the emission spectra of D2 obtained during our 
studies of translational and rotational temperatures in hydrogen and deuterium containing 
plasmas {3^. Detailed description of the experimental setup was reported elsewhere 



Capillary arc discharge lamps DDS-30 described in [48| have been used as a light source. 
They were filled with about 6 Torr of the spectroscopically pure D2 + H2 (9:1) mixture. 
The range of the discharge current was from 50 to 300 mA (current densities j = 1.6 4- 10 
A/cm^). The light from the axis of plasma inside the capillary was directly focused by an 
achromatic lens on the entrance slit of the Czerny- Turner type 1 m double monochromator 
(Jobin Yvon, UIOOO). The intensity distribution in the focal plane of the spectrometer was 
recorded by cooled CCD matrix detector of the Optical Multichannel Analyser (Princeton 
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Appl. Res., OMA- Vision-CCD System). 

For identification of tlie D2 spectral lines in our spectrum we used tlie assignments and 
wavelengtli values from IGj. Tfiose values show certain spread around monotonic depen- 
dence of the wavenumber from the distance along the direction of dispersion in the focal 
plane of the spectrograph. The dispersion function of the spectrograph was obtained by 
the polynomial least squares fitting of the wavelengths versus the distance for strong un- 
blended lines. Such wavelength calibration allow us to get new experimental values for the 
wavenumbers of 125 questionable spectral lines. The new values differs from those reported 
in not only because of the smoothing procedure. In contrast to previous works we used 
digital intensity registration providing linear response of the CCD detector. That gave us 
the opportunity of digital deconvolution of the recorded spectra to resolve the major part 
of blended lines. The observed line profiles were determined mainly by Doppler and instru- 
mental broadening. For strong unblended lines they were close to Gaussian shape except of 
insignificant far wings (see e. g. bottom graph in figure [5l). Therefore, the parts of spectra 
in the neighbourhood of the blended questionable lines under the study were approximated 
by the superpositions of certain number of lines having Gaussian profiles with fixed half 
width and adjustable intensity and wavelength values. Then we obtained new values for the 
maximum intensities of D2 lines. Typical example is shown in the figure O One may see 
that the new values of the rovibronic transition wavenumbers are a little bit different from 
those reported in 16|. In our experimental conditions the intensities of 10 questionable lines 
were too weak for detecting and unambiguous separation by the deconvolution. 

The uncertainty of our experimental data is mainly determined by errors of the decon- 
volution process. It should be underlined, that most of the questionable lines are blended 
lines of relatively low intensity. In that case it is not easy to obtain reliable RMS estimates 
for experimental errors. Therefore we used statistical method for derivation of experimental 
uncertainty described above. The new wavenumber values were included into the input data 
set obtained in the last stage of previous section as a separate (20-th) sample of uniformly 
precise experimental data. After several iterations the root mean square estimate of exper- 
imental uncertainty 0.06 cm^^ was obtained without any evidence for existence of outliers. 
The experimental $(0 function for the sample of new wavenumbers for questionable lines 
is shown in Figure [6] together with corresponding frequency diagram. One may see that 
the distribution function of random errors is close to the normal distribution. Moreover 
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new wavenumber values do not contradict to the rest of the data set and may be used for 
determination of optimal values of rovibronic levels. 

The new wavenumber values of questionable lines obtained in the present work are listed 
in table [T] together with those from [l^ . One may see that in most cases the differences 
between new and old data are significant. Most often small differences (< 0.05 cm~^) are 
probably caused by errors of reading from photographic plates and/or by round-up errors, 
when the vacuum wavenumbers were calculated from the wavelengths measured in air. In our 
case they were eliminated by digital reading of the intensity profiles and smoothing procedure 
of the determination of the dispersion curve. The differences higher than 0.05 cm~^ are 
caused by blending (by shifts of the intensity maxima or by disappearance of a weak line 
within the profile of the strong one, see e.g. figure [5]), misprints and wrong assignments. The 
differences (0-C) between Observed wavenumber values and those Calculated as differences 
of corresponding optimal energy level values obtained in the present work are also shown in 
table [H One may see that in the framework of Rydberg-Ritz principle the new experimental 
data for 125 questionable lines are in much better co-ordination with the 3588 wavenumbers 
of other lines than those reported in . 



V. RESULTS AND DISCUSSION 



Information on the amount of experimental data reported in original papers and the 
amount of outliers revealed by statistical analysis is presented in table [Tll 

After all, the input data set consists of 3713 wavenumber values (3588 old and 125 new) 
divided into 20 subsets of uniformly precise data with known RMS estimates of experimental 
uncertainties obtained in the present work by pure statistical analysis, i.e. independent 
from the estimates reported in original papers. This data set was used for determination 
of optimal values of all studied experimentally triplet rovibronic levels by minimizing (I3l). 
The empirical function $(0 thus obtained is presented in figure [3 One may see that the 
final distribution function -F(^) is close to the normal distribution -fo(0- Small deviations 
from the linear plot are caused by insufficient statistics for the data showing too high values 
of the weighted deviations (without any contradiction with the Chauvenet's criterion). 
Thus, optimal values of the 1050 rovibronic energy levels have been obtained from the 3713 
values of the experimental wavenumbers. All the energy levels were obtained relative to the 
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lowest vibro-rotational level (f = 0, = 0) of the electronic state. 

The shift (33.631 ± 0.004) cm~^ between uncoupled levels of ortho- and paradeuterium 
was obtained by the least squares analysis of the f = 0, = -i- 18 levels with 

odd and even values of the rotational quantum number A^. According to [l^ the difference 
Eaoi — Eaoo = (33.62 ± 0.10) cm~^. Both values coincide within error bars, but our value is 
obtained with much higher precision. Taking into account the values of the SD uncertainties 
of our optimal level values (see table HTB . this increase of precision in obtaining the shift 
value is significant. 

Optimal values of the rovibronic levels for the triplet electronic states corresponding to 
principle quantum numbers = 2, 3 of the united atom are listed in table lllli Radiative 
transitions between those levels are located in visible part of spectra, which is most conve- 



he expanded 



46| . For each 



nient for spectroscopic studies of non-equilibrium gases and plasmas 
version of the table UTTl including all 1050 optimal level values is presented in 
rovibronic level we introduced into the table IIIII the standard deviation of the empirical 
determination (in brackets), the number of spectral lines n^, which originate or terminate 
on that level, and the difference AE between energy level values reported in 161] and those 
obtained in the present work. 

Several interesting things could be seen from the table IIIII namely: 

i) In many cases the number of the lines, directly connected with certain level, is high 
enough for statistics. In some cases there are very few lines. In that cases the levels are 
obtained with less accuracy. It is important to stress, that the SD uncertainties of empirical 
determination of optimal level values partly reflect the influence of rii, on the SD values. 
This is the trend (not straightforward dependence) because there are more factors of the 
influence. The levels with n^, < 4 should be considered as unreliable. The SD of such levels 
are mainly determined by experimental uncertainties of few lines directly connected with 
those levels. Additional experimental studies are needed to obtain more reliable results. 

ii) The differences AE are generally less, than the value 0.05 cm~^, reported in [l^ 
as SD uncertainty of the energy levels. This confirms our assumption that the unique 



16| is only an upper-limit estimate of real 



SD uncertainty for all the levels declared in 
uncertainties of various levels. 

iii) The vast majority of the AE are much higher than SD error bars of our optimal 
values, which normally are within the range of 0.004 0.03 cm~^ depending on the value of 
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rotational and vibrational quantum numbers. Thus, the deviations of the data reported in 

n 

|16l | from the optimal level values obtained in the present work are significant. 

Currently it is not possible to find out and list reasons for those deviations for all and every 
levels, because IGj does not contain information about details of calculating the common 
differences in various bands. We may presume that they appeared as a result of following 
reasons: errors of reading from photo plates, round up errors in calculating the wavenumbers 
from measured wavelengths, lack of deconvolution of blended lines, subjective preferences 
in eliminating some of controversial experimental data, some additional assumptions (like 
neglecting A — doubling for some states), transfer ^d accumulation of experimental errors 
in the method of common differences adopted in [16] for wavenumber analysis of separate 
bands, quite subjective multi-step procedure of the analysis of the band systems. It is clear 
that for various rovibronic states all these factors may multiply or compensate each other. 

Sometimes the nature of the differences AE may be easily recognized. One typical ex- 
ample is illustrated in figure [HI in which the differences AE are shown for rotational levels 
in the ground vibrational states {v = 0) of the (d), c^H^ (c), and j'^A^ (b) electronic 

states. The error bars represent our estimates (one SD) of the empirical determination of the 
optimal energy level values obtained by pure statistical approach in the present work. One 
may see that, as it was already mentioned, quite often the deviations AE are less than 0.05 
cm~^, but significantly higher than SD uncertainties of optimal level values (figure [8]^d)). In 
figure [H](b) and figure [H](c) one may see interesting behaviour of the deviation AE — peri- 
odical changes for rotational level with odd and even values of (for uncoupled levels of 
ortho- and para- molecules) . Similar effect is observed for many other electronic- vibrational 
states including the case shown in figure [8](d). 

Small alternations of the deviations AE (less than 0.1 cm~^) are connected with the 



difference in the values o 
the present work and in 



the shift between levels of orto- and para- molecules obtained in 



IGj. More pronounced cases of the alternations appear as a result 
of experimental errors due to non-optimal method of energy level derivation used in [l^ — 
the sequential, multi-step procedure based on calculating common differences. The error 
appeared in one common difference due to one spurious experimental wavenumber value is 
transferred from previous energy level value to the next one. The mechanism is valid also 
for small random errors with various final results, because subsequent pairs of lines within 
a band may have experimental errors of various values and sings. Sometimes random errors 
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may partially compensate each other, but in principle this method has the possibility of an 
accumulation of experimental wavenumber errors and transferring them into the errors of 
the empirical values of rovibronic energy levels. 

The typical example of the alternation appearance caused by spurious experimental 
wavenumber values is illustrated in figure [^a,b,c). They demonstrate the influence of two 
s,gn,flca„t errors m wave.umbers of the fl-branch li,res of the (0-0) band of the fA,-^U, 
electronic transition on the empirical values of the rotational levels obtained in |16|] for the 
upper , V = and lower c^H", electronic- vibrational states of the band. 

n 

Excellent agreement between the observed and calculated wavenumbers reported in [IS] 
seen in figure [Hl(a) is the evidence that the wavenumbers of the R4 and -R12 spectral lines 
were used in determination of the level values. Experimental wavenumbers reported in 



ts of our 



16|. The 



16l | for both lines were classified as outliers by our statistical analysis. Resu 
independent experiments confirmed that this is caused by occasional errors in 
wavenumber values from [l^ for the lines RA and -R12 are underestimated and overestimated, 
respectively (see figure [H]( a) and table [T]). From the figure [H^b,c) one may see how those 
experimental errors are transferred into the wrong values of the empirical energy levels via 
sequential adding of the common differences A2(A^) and A2(A^) calculated for upper and 
lower rovibronic levels. The cases of appearance, accumulation and partial compensation of 
experimental errors may be found in figure [8]J^b,c). 

Takin g in to account the results of our previous studies of some triplet states of H2 
molecule 26|| it is easy to expect that the Dunham series expansions with empirical molecu- 



lar constants from 



15, 



251 ] have rather limited ability to describe rovibronic levels of the D2 



isotopomer. These expectations were confirmed by the results of the present work. Typical 
example of the contradictions is illustrated in figure [9] representing the differences AE be- 
tween optimal values of the energy levels obtained in the present work and those calculated 



with molecular constants from 



15|, 



251 ] for various rotational and vibrational levels of the 



a^S+ electronic state of D2 molecule. One may see that the deviations show strong de- 
pendencies on vibrational and rotational quantum numbers, achieving values of about 100 
cm~^. It should be stressed that according to our previous studies of the potential curves [s^ 
the electronic state is almost free from non-adiabatic effect of perturbations. Higher 

electronic states are mainly perturbed due to electronic-rotational and electronic-vibrational 
interactions, and simple Dunham expansions are non applicable. 



For determination of precise absolute values of the triplet rovibronic levels of D2 it is 
desirable to make pure statistical analysis of all currently available experimental data on 
the wavenumbers of singlet rovibronic transitions and to use the data of anticrossing spec- 



troscopy 



49, 



5G, 



5l| for establishing the link between singlet and triplet state. The approach 



used in the present work is most appropriate for solving this problem in the future. Nowa- 
days the absolute values of the triplet rovibronic levels (relative to the X^Eg+,f = OA^ = 
ground rovibronic state) may be obtained only by adding an absolute value -Eaoo of energy 
of the fl^S^, f = = rovibronic level to the relative energy level values obtained in the 
present work. According to 16] ^^00 = 95348.2(4) cm'^ while in E^oo = 95348.3(1) 
cm~^. These values coincide within the uncertainties reported in original papers and may 
be used. The precision of such absolute calibration should be better than 1 cm~^. 

Our relative values of triplet rovibronic levels of D2 molecule may be used not only in 
comparisons with results of non-empirical calculations. They are recommended for accurate 
calculating the wavenumbers of triplet-triplet optical transitions for computer simulations 
of emission and absorption spectra in applied spectroscopy of non-equilibrium gases and 
plasmas. Moreover, the high precision of optimal level values achieved in the present work 
provides the opportunity to expand existing identification of triplet rovibronic lines by accu- 
rate calculating wavenumbers of all currently unassigned rovibronic transitions allowed by 
the selection rules for electric dipole transitions and detecting corresponding lines in exper- 
imental spectra. From the Grotrian diagram of the D2 molecule it may be seen that several 
new band systems may be discovered and investigated. 
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FIG. 1: Empirical function $(^) obtained when all published experimental data with RMS uncer- 
tainties equal to the error estimates reported in original papers. The onset represents the central 
part of the in larger scale. 
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FIG. 2: Empirical functions $*^"*)(^) obtained in different m-th iterations (a, b, c, and d correspond 
to m = 1, 2, 3, and 4) of the statistical analysis of the sample of experimental data representing 
wavenumbers of R- and P- branch lines of the (0 — 1) band of the e^E+ — a^S+ electronic transition 
of D2. 
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FIG. 3: The amount of the wavenumber values rii, within various enlarged samples of the uniformly 
precise experimental data. 
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FIG. 4: Empirical functions $(^) representing the error distributions within the first six enlarged 
samples corresponding to the intervals (Oh-0.01), (0.01 H- 0.02), (0.02H-0.03), (0.03H-0.04), (0.04-=- 
0.05), (0.05 0.06) (from left to right, from top to bottom). The values of the experimental 

uncertainty obtained in the present work for every sample are shown near appropriate graphs 
together with SD of their statistical determination (in brackets). 
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FIG. 5: Fragments of the D2 spectrum in the neighbourhood of the questionable lines. Experimen- 
tal intensity J in relative units is shown by open circles. Dotted lines represent Gaussian profiles 
for separate lines obtained by deconvolution, while solid line corresponds to the total inten sity 



16| 



obtained by summing over the components. The wavenumber values from the Appendix C of 
are shown and marked by dash-dot vertical lines. Pieces of bold vertical lines denote maxima of 
the separate line profiles. Bold underlined wavenumbers represent the new wavenumber values for 
the questionable lines obtained in the present work. 
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FIG. 6: Empirical function $(^) for the sample containing new wavenumber values for questionable 
lines. The frequency diagram N versus ^ is constructed with uniform width intervals equal to 0.4. 
N represents the number of the deviations ^ij within each interval. The columns show experimental 
results while the solid line with circles corresponds to the normal distribution with zero mean and 
variance equal to unity. 
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FIG. 7: Empirical function $(^) obtained for the final set of all available experimental data (see 
the text) with RMS errors estimates obtained in the present work by the statistical analysis. 
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FIG. 8: The differences (O — C) between observed wavenumbers from [2] and those calculated 
as the differences of empirical levels from jl6] (1) and of the optimal level values (2) for i?-branch 
lines of the (0 — 0) band of the J^A~ — c^n~ electronic transition (a). The error bars in figure 
(a) represent the sum of standard deviations for optimal values of upper and lower levels and for 

n 

the wavenumber value (obtained in the present work for the set of experimental data from [la ] 
belonging to this band). The differences AE between the level values from {igI and the optimal 
level values (Table lin|) for various rotational levels of the j^A~, v = (b); c^H", v = (c); and 
a^S^, V = (d) electronic- vibrational states of D2 molecule. Spectral lines of the R- and P- 
branches of the (0-0) band of the J'^ A~ — c^n~ electronic transition connecting some of the levels 
are shown by vertical lines (dashed for i2-branch and solid for P-branch) with standard labels. 
The error bars in figures (b), (c) and (d) represent SD of the optimal energy levels obtained in the 
present work. Horizontal solid lines represent zero level, and horizontal dash-dotted lines — error 
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FIG. 9: The differences AE between optimal values of the energy levels obtained in the present 
work and those calculated with molecular constants from for various rotational and vibrational 
levels of the electronic state of D2 molecule. 
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TABLE I: Experimental values of the wavenumbers (in cm^^) of some spectral lines of D2 obtained 
in the present work (P.W.)and those reported in IgI. 0-C denotes the differences between Observed 
wavenumber values and those Calculated as differences of corresponding optimal energy level values 
obtained in the present work. 

Electronic transition Band Line [16] 0-C P.W. 0-C 

e^j:+ - a3S+ 2-0 Pll 13751.08 -0.14 13751.29 0.07 

3- R6 15868.98 -0.05 15869.01 -0.02 

R9 15705.21 0.23 15704.98 0.00 

4- 1 P8 14862.89 -0.08 14862.95 -0.02 

5- 2 P6 14524.31 0.07 14524.26 0.02 

6- 2 R4 15977.48 -0.06 15977.56 0.02 

6- 3 P9 13703.95 0.09 13703.87 0.01 

7- 3 P5 15107.67 -0.04 15107.71 0.00 

8- 3 RO 16354.21 0.05 16354.19 0.03 

8- 4 P5 14468.77 -0.10 14468.86 -0.01 

R3 14727.33 -0.06 14727.38 -0.01 

9- 5 P4 13861.27 -0.19 13861.45 -0.01 

R4 14002.82 -0.11 14002.91 -0.02 

10- 4 P4 16124.14 0.25 16123.94 0.05 
- a^T.+ 1-1 R2 20346.49 0.07 20346.42 0.00 

1- 2 RO 18564.54 0.21 18564.19 -0.14 

2- 1 P3 21653.13 0.20 21653.04 0.11 
2-3 P3 18242.42 0.09 18242.40 0.07 

d^Uu-a^^j 0-1 R9 15019.54 0.10 15019.44 0.00 

1-0 P2 18207.70 0.10 18207.60 0.00 

P5 18065.38 0.06 18065.36 0.04 

1-1 Ql 16460.87 -0.07 16460.93 -0.01 

QIO 16267.60 -0.05 16267.64 -0.01 

1- 2 QIO 14582.86 0.03 14582.84 0.01 

2- Ql 19823.53 0.09 19823.45 0.01 
2-1 Rl 18065.70 0.07 18065.68 0.05 

R5 18109.66 0.04 18109.64 0.02 

R6 18107.89 -0.04 18107.97 0.04 
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table IB (Continued.) 
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TABLE II: The amounts of various kinds of the experimental data reported in various original 
papers and obtained in the present work (P.W.): total reported, useable for the statistical analysis, 
excluded as the outliers, used for determination of optimal level values, and the rate of used lines 

in %. 

Reference [16|[32][3l][28][34][33l[29l[35] P-W. Total 
reported 3117 350 285 81 37 83 3 1 125 4082 
useable 3053 350 285 81 37 12 3 1 125 3947 
excluded 179 18 18 5 4 10 234 
used 2874 332 267 76 33 2 3 1 125 3713 
% 94.1 94.9 93.7 93.8 89.2 16.7 100 100 100 94.1 
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TABLE III: Optimal values EnvN of rovibronic energy levels (in cm^^) for various triplet elec- 
tronic states of D2 molecule obtained in the present work. The uncertainties of the EnvN value 
determination (one SD) are shown in brackets in units of last significant digit, riy — the number 
of various spectral lines directly connected with certain level. A.E is the difference between energy 
level values obtained in the present work and those reported in 





V — 







V — 


1 






V = 


2 




N 


EnvN 




A F, 


EnvN 






EnvN 


/ iy 




n 

u 


n 

u 


91 


n 
u 


1 k1 A 


9"^ 
zo 


00^ 
-u.uuo 


oOOU 


1 ZD(D j 


9f^ 
zu 


01 R 

-U.UlU 


1 

1 


OO.UOl y± J 


'4:0 


-U.Ul 1 


1 RA7 '^(\A(R\ 


00 


01 zl 


3592 


036(8j 


ou 


01 
-u.uiu 


9 


1 no 77n((\^ 


oz 


01 
-u.uxu 


1 Ql 9 1 791'^^ 


uu 


01 9 


oD04 


( 6[o ) 


70 
/ u 


00"^ 
-u.uuo 


Q 


901 inoi'Q^ 
zui. luy i^y ) 




OOQ 


900Q 0A9('^'\ 


Do 


009 

-u.uuz 


( 4o 


nQ/i ( Q\ 


(\A 


00/1 
-u.uU^ 


A 

4: 






01 
-u.uxo 




^7 


00^ 
-u.uuo 


00 ( / 




(\A 
u^ 


00^ 
-u.uuo 




4QQ R^f)( 1 Cl\ 




090 


99Q7 '^^^^('Q'i 


■4:1^ 


00^ 
-u.uuo 


4UZD 


z4o(y j 




00"^ 
-u.uuo 


a 
u 


(\Q7 nQ9l'7'\ 
1 .uyz ^ / j 




019 


9zlR7 7/lki'fi'\ 


49 

'4:Z 


009 
u.uuz 


4zU9 


894( / j 


AA 


OOfi 
u.uuu 


7 


09"^ 9fini' 1 1 


1 Q 


090 


9707 Q(\7( 1 0^1 


97 

Z 1 


007 
-u.uu / 


4422 


319(10) 


9"^ 
zo 


001 

U.UUl 


Q 



iioo.O 1 oyo j 


1 Q 


01 


9Q'^7 990(^7^ 
zyo < .zzui^ i j 


97 
Z 1 


090 

u.uzu 


4662 


756(7) 


99 
zz 


01 /I 


9 


1471.063(13) 


11 


-0.043 


3234.634(11) 


21 


-0.004 


4930 


308(11) 


14 


-0.008 


10 


1786.713(9) 


12 


-0.023 


3539.206(8) 


15 


0.024 


5224 


035(9) 


10 


0.015 


11 


2129.524(17) 


6 


-0.054 


3869.902(14) 


6 


0.008 


5542 


94(2) 


4 


-0.010 


12 


2498.252(13) 


6 


-0.042 


4225.651(16) 


3 


0.069 


5885 


90(3) 


3 


0.010 


13 


2891.75(2) 


5 


-0.019 


4605.20(2) 


3 


-0.010 


6251 


85(3) 


2 


0.040 


14 


3308.86(2) 


5 


-0.050 


5007.43(3) 


3 


-0.050 


6639 


59(3) 


2 


0.020 


15 


3748.44(3) 


5 


-0.010 


5431.37(3) 


3 


0.010 


7048 


04(3) 


2 


0.030 


16 


4208.73(3) 


5 


0.070 


5875.16(4) 


2 


0.050 


7475 


59(17) 


2 


0.050 


17 


4689.10(3) 


4 


0.010 


6338.32(4) 


2 


0.040 










18 


5187.38(4) 


3 


0.060 


6818.59(4) 


2 


0.080 











v = 3 v = 4 v = 5 

N EnvN n-i, AE EnvN ny AE EnvN ny AE 

5238.994(6) 27 -0.014 6851.077(6) 21 -0.007 8398.070(7) 14 -0.000 

1 5269.198(8) 54 -0.008 6880.190(9) 48 0.000 8426.100(9) 37 0.000 

2 5329.503(6) 64 -0.013 6938.313(6) 57 -0.013 8482.065(6) 37 0.015 



table [ml (Continued.) 





V = 6 


V = 7 


V = i 






N EnvN 


/S.E EnvN IT'U 


AE EnvN 




AE 


9880.806(11) 10 


-0.016 11299.845(16) 6 


-0.045 12655.467(18) 


5 


-0.007 


1 9907.775(10) 26 


0.005 11325.784(11) 18 


-0.004 12680.356(12) 12 


-0.036 


2 9961.617(7) 32 


-0.017 11377.500(8) 18 


0.000 12729.972(9) 


13 


-0.052 


3 10042.093(10) 28 


-0.013 11454.821(12) 14 


-0.010 12804.129(13) 14 


-0.059 


4 10148.876(7) 28 


-0.026 11557.399(9) 14 


n nnQ i9Qn9 ^d'^i'q'i 

u.uuy ijLicjyj^ .oyjoyij J 




-0.143 


5 10281.509(10) 25 


-0.069 11684.766(13) 8 


-0.036 13024.059(14) 


4 


-0.049 


6 10439.463(8) 20 


-0.033 11836.493(10) 7 


-0 033 13170 42(5) 


2 


-0.000 


7 10622.101(13) 12 


-0.071 12011.941(19) 3 


-0.051 






8 10828.794(13) 7 


-0.084 12210.76(3) 2 


-0.060 






9 11057.80(2) 2 


0.010 








10 










11 11594.7(3) 2 










12 11894.2(3) 2 










V = 9 


V = 10 








N EnvN nu 


AE EnvN riy 


AE 






13947.39(4) 1 


0.000 








1 13971.355(15) 4 


-0.035 15198.29(5) 2 


-O.OoO 






2 14018.883(12) 4 


-0.013 15243.76(5) 2 


-0.020 






3 14089.876(15) 3 


0.024 15311.62(5) 2 


-0.040 






4 14184.038(13) 4 


-0.088 15401.72(7) 1 


0.010 






5 14300.933(17) 3 


0.017 15513.75(5) 2 


-0.050 






6 14440.112(13) 1 


0.028 


















V = 


V = 1 


V = : 


2 




N EnvN 


AE EnvN 


AE EnvN 




AE 


1 -132.611(18) 12 


0.021 1548.794(14) 9 


0.036 3168.223(11) 


8 


0.087 


2 -71.553(12) 17 


-0.007 1607.788(9) 15 


0.032 3225.211(8) 


12 


0.079 



45 



table [ml (Continued.) 







V = 3 


V = 4 


V = 5 


N 


EnvN n-jy AE 


EnvN ny AE 


EnvN ny AE 



1 4727.022(9) 7 0.108 6226.366(19) 3 0.144 7667.225(13) 3 0.125 

2 4782.112(9) 10 0.038 6279.454(11) 7 0.096 7718.429(10) 2 0.141 

3 4864.307(13) 11 0.073 6358.796(14) 4 0.084 7794.944(13) 2 0.116 

4 4973.456(10) 12 0.054 6464.117(13) 4 0.063 7896.475(10) 2 -0.035 

5 5109.097(16) 9 0.093 6594.990(18) 4 0.050 

6 5270.732(15) 5 0.028 

7 5457.84(2) 5 -0.010 

8 5669.46(2) 4 0.030 

9 5905.35(3) 1 -0.020 
10 6164.52(3) 1 0.080 



V = 6 



N EnvN rLy 


AE 








1 9050.44(2) 1 










2 9099.750(17) 1 










3 9173.43(2) 1 










4 9271.249(18) 1 




















V = 


V = 1 


V = 


2 




N EnvN n-u 


AE EnvN riy 


AE EnvN 




AE 


1 -132.586(12) 22 


-0.004 1548.813(9) 16 


0.017 3168.273(7) 


13 


0.037 


2 -71.584(14) 26 


0.014 1607.790(12) 19 


0.040 3225.247(10) 17 


0.123 


3 19.662(11) 28 


-0.032 1695.990(8) 21 


0.040 3310.490(9) 


16 


0.030 


4 140.727(15) 27 


-0.026 1813.023(13) 18 


0.047 3423.541(11) 


14 


0.079 


5 291.217(12) 24 


-0.087 1958.482(10) 18 


0.038 3564.085(8) 


14 


0.095 


6 470.542(16) 19 


-0.022 2131.814(15) 14 


0.016 3731.563(12) 


8 


0.057 


7 678.114(15) 15 


-0.094 2332.393(14) 9 ■ 


-0.003 3925.279(18) 


9 


0.051 


8 913.12(2) 10 


0.022 2559.53(2) 5 


0.060 4144.65(2) 


4 


0.140 



46 

table mi (Continued.) 



{lsa2pTT)c^U: 



V = 3 






V = 4 




V = 5 






N EnvN 






EnvN ny 


/S.E EnvN fly 


AE 


1 4727.078(7) 


14 


0.052 


6226.429(10) 


7 


0.081 7667.297(13) 


2 


0.053 


2 4782.081(12) 16 


0.109 


6279.553(10) 


9 


0.007 7718.544(13) 


2 


0.016 


3 4864.446(8) 


15 


0.014 


6359.007(11) 


7 


0.133 7795.177(10) 


1 




4 4973.606(14) 


9 


0.084 


6464 401 I'M"! 


8 


079 7896 833(13") 


1 


- 


5 5109.368(10) 10 


0.092 


6595.379(13) 


6 


-0 009 8023 166fl0) 


1 


- 


6 5271.03(2) 


5 


0.110 


6751.361(18) 


4 


0.059 






7 5458.145(18) 


7 


0.035 












8 5670.02(3) 


5 


-0.030 












9 5905.61(2) 


3 


0.070 












V = 6 




V = 7 










N EnvN 


riu 


AE 


EnvN ny 


AE 






1 9050.47(2) 


1 


- 


10376.62(2) 


1 








2 9099.88(2) 


1 




10424.20(2) 


1 








3 9173.710(17) 


1 




10495.262(18) 


1 








4 9271.63(2) 


1 




10589.37(2) 


1 








5 9393.409(18) 


1 




10706.613(18) 


1 














(lscj3scr)/i3E+ 








V = 2 






V = 3 




V = 4 






N EnvN 


riy 


AE 


EnvN ny 


AE EnvN fly 


AE 


20071.893(17) 


1 


0.037 21588.238(10) 


1 


0.052 23042.197(13) 


2 


0.083 


1 20101.320(16) 


2 


0.100 21615.731(12) 


2 


0.099 23068.46(2) 


2 


0.143 


2 20160.107(13) 


2 


0.053 21670.798(8) 


3 


0.032 23120.982(10) 


4 


0.088 


3 20247.847(15) 


3 


0.113 21752.929(13) 


3 


0.061 23199.357(13) 


4 


0.083 


4 20364.077(13) 


2 


0.063 21861.839(9) 


3 


0.051 23303.289(13) 


3 


-0.009 


5 










23432.221(16) 


3 


0.059 
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table [ml (Continued.) 

(lscr3p(T)e^S+ 



v=0 v=l v=2 



N EnvN 




AE 


EnvN 




AE 


EnvN 




AE 


11649.190(14) 


5 


-0.010 


13137.502(10) 


6 


-0.012 


14558.258(12) 


6 


0.002 


1 11676.174(8) 


12 


-0.004 


13163.431(6) 


12 


-0.011 


14583.137(7) 


14 


0.003 


2 11730.037(13) 12 


-0.007 


13215.170(9) 


15 


-0.010 


14632.770(11) 


13 


0.000 


3 11810.600(9) 


17 


-0.020 


13292.570(7) 


16 


-0.010 


14707.003(8) 


15 


-0.003 


4 11917.524(12) 12 


-0.004 


13395.293(10) 


11 


-0.003 


14805.515(11) 


14 


0.015 


n 1 20.^0 48.5 ('q') 


12 


0.005 


13523.013(7) 


12 


0.007 


14927.986(8) 


14 


-0.006 


6 1 2208 941 (^ 3"! 


s 


-0.021 


13675.191(10) 


12 


-0.001 


15073.923(11) 


12 


-0.013 


7 12392.375(9) 


8 


0.005 


13851.369(8) 


11 


0.011 


15242.814(8) 


13 


0.016 


8 12600.117(13) 


6 


-0.017 


14050.868(11) 


9 


-0.018 


15434.079(13) 


8 


-0.039 







-0.030 


14272.987(8) 


10 


0.003 


15646.946(10) 


8 


0.024 


-LU J.OuO(J.UiJO^^lU } 


c; 
(J 


-0.043 


14517.021(13) 


6 


-0.021 


15880.745(15) 


5 


0.005 


11 13361.853(12) 


4 


-0.013 


14782.091(10) 


5 


-0.001 


16134.676(15) 


3 


0.044 


12 13659.139(19) 


4 


-0.019 


15067.365(16) 


3 


0.015 


16407.97(10) 


2 


-0.060 


13 13976.603(17) 


3 


-0.043 


15371.957(16) 


5 


-0.007 


16696.56(5) 


2 


-0.030 


14 14313.27(2) 


3 


0.000 


15694.80(2) 


4 


0.030 


17003.73(4) 


2 


0.010 


15 14668.18(3) 


2 


-0.050 


16035.06(3) 


4 


0.050 


17326.48(5) 


2 


0.070 


lb 15040.49(3) 


2 


0.010 


16391.90(3) 


4 


0.010 


17675.65(3) 


1 


0.010 


17 15428.58(4) 


2 


0.050 


16763.73(4) 


4 


0.060 








18 15832.13(3) 


2 


0.020 


17150.23(4) 


2 


0.040 








19 16249.09(4) 


2 


0.060 


17549.32(5) 


1 


0.080 








V = c 


5 




V = 4 


I 




V = I 






N EnvN 


riu 


AE 


EnvN 


riy 


AE 


EnvN 






15911.178(10) 


7 


-0.018 17195.132(10) 10 


-0.002 18408.097(10) 


7 


-0.007 


1 15934.967(6) 


19 


-0.057 17217.850(6) 


17 


0.000 18429.790(6) 


16 


-0.010 


2 15982.471(9) 


20 


-0.001 17263.194(8) 


17 


-0.004 18473.084(9) 


15 


-0.006 


3 16053.509(6) 


19 


-0.009 17331.018(6) 


19 


0.002 18537.934(7) 


14 


-0.034 


4 16147.768(9) 


19 


0.012 17420.989(9) 


18 


-0.008 18624.360(9) 


14 


0.010 


5 16264.946(6) 


18 


-0.006 17532.827(6) 


17 


-0.007 18734.130(6) 


14 


-0.010 



table [ml (Continued.) 



(lscr3po")e^S+ 


V = 6 


V = 7 


V = i 


i 




N EnvN 


/S.E EnvN f^u 


AE EnvN 




AE 


19546.641(11) 7 


-0.001 20605.338(10) 8 


-0.008 21575.310(13) 


9 


-0.030 


1 19567.047(7) 16 


-0.007 20624.515(6) 15 


-0.015 21593.153(7) 


19 


-0.013 


2 19607.728(10) 16 


-0.008 20662.787(9) 16 


-0.017 21628.745(10) 19 


-0.005 


3 19668.569(7) 16 


0.001 20720.289(6) 14 


-0 009 21681 957l'7') 


20 


-0.007 


4 19749.319(10) 16 


-0.019 20795.613(9) 12 


-0 023 21752 568fll) 


18 


-0.038 


5 19849.642(8) 15 


-0.012 20889.773(6) 14 


-0.003 21840.510(8) 


16 


-0.010 


6 19969.098(11) 14 


-0.028 21001.748(10) 9 


-0 008 21946 349fl2) 


7 


-0.029 


7 20107.185(7) 12 


0.005 21131.071(7) 9 


-0 001 22074 088fl4) 


2 


0.002 


8 20263.354(13) 8 


0.006 








9 


21440.901(17) 1 


0.029 






V = 9 


V = 10 








N EnvN f^v 


/S.E EnvN 








22442.451(11) 8 


-0.011 23181.37(2) 4 


n nnzL 






1 22458.482(7) 17 


-0.002 23195.302(13) 11 


-0 099 






2 22490.446(10) 16 


-0.006 23223.064(16) 10 


-0.024 






3 22538.192(7) 18 


-0.002 23264.357(13) 11 


-0.017 






4 22601.388(10) 18 


-0.038 23318.792(19) 6 


-0.042 






5 22679.660(7) 14 


-0.040 23385.89(2) 3 


-0.030 






6 22772.422(11) 10 


-0.022 








7 22879.133(11) 3 


0.037 


















V — n 

V — u 


V L 


V = ^ 






N EnvN fly 


AE EnvN flu 


AE EnvN 


riu 


AE 


1 16695.750(14) 4 


-0.010 18308.303(11) 7 


-0.033 19857.067(10) 


8 


0.003 


2 16755.305(11) 5 


-0.005 18365.680(8) 5 


-0.010 19912.363(7) 


8 


-0.023 


3 16844.286(15) 4 


-0.006 18451.471(11) 7 


-0.001 19995.024(10) 


8 


-0.014 


4 16962.446(11) 4 


-0.006 18565.354(8) 7 


-0.004 20104.741(7) 


8 


-0.011 



49 

table [ml (Continued.) 



V = 3 




V = 4 




V = f 






N EnvN riy 


Ai? EnvN 


AE EnvN 






1 21343.645(12) 


9 


0.005 22769.442(12) 


8 


-0.012 24135.655(10) 10 


0.005 


2 21396.902(9) 


8 


-0.002 22820.678(8) 


8 


002 941 84 948l'7) 


1 


-0.008 


3 21476.477(12) 


8 


-0.007 22897.296(12) 


7 


-0 006 24258 610n0) 


9 


0.000 


4 21582.138(9) 


8 


-0.008 22998.963(8) 


8 


-0 013 24356 M7(7) 


9 


-0.017 


5 21713.423(13) 


7 


-0.013 23125.264(12) 


8 


-0 004 24477 809(1 1 ) 


7 


-0.029 


6 21869.803(9) 


8 


0.057 23275.747(9) 


7 


0.003 24622.462(9) 


9 


-0.012 


7 22050.667(15) 


4 


-0.007 23449.760(17) 


5 


-0.020 24789.747(12) 


7 


-0.007 


8 22255.379(12) 


4 


0.011 23646.709(11) 


6 




D 


-0.010 


Q 22483 185 f 19) 


2 


-0 035 23865 2(^(2^ 


2 




o 


0.010 


10 22733.258(15) 


3 


0.002 24106.31(3) 


2 


0.020 25420.32(3) 


2 


0.040 


11 23004.67(5) 


2 


-0.050 24367.28(6) 


1 


-0.030 






12 23293.59(6) 


2 


-0.010 24645.28(7) 


1 


-0.010 






13 23603.42(6) 


1 


0.030 










14 23930.61(6) 


1 


0.090 










V = 6 




V = 7 




V = ^ 






N EnvN riy 


/S.E EnvN 1 




AE EnvN 




AE 


1 25443.364(10) 


9 


-0.004 26693.448(12) 


6 


-0.008 27886.595(13) 


6 


-0.005 


2 25490.706(7) 


9 


-0.006 26738.891(9) 


6 


-u.uzi z I you. lozi^y j 


D 


0.008 


3 25561.485(11) 


8 


-0 005 26806 814('12') 


6 




u 


0.035 


4 25655.396(7) 


8 


-0 006 26896 934(9) 


5 


-0.054 28081.598(10) 


6 


-0.028 


5 25772.049(11) 


9 


n ni Q 97nns 
-u.uiy z(UUo.oo(^zj 


c; 



-0.026 28188.852(15) 


6 


-0.072 


6 25910.990(8) 


7 


-0.020 27142.454(10) 


3 


-0.024 28316.512(10) 


5 


0.008 


7 26071.762(19) 


3 


-0.082 27296.310(15) 


2 


-0.070 






8 26253.93(3) 


2 


-0.080 27470.81(4) 


1 


-0.060 






V = 9 




V = 10 










N EnvN ny 


AE EnvN f 




AE 






1 29023.33(5) 


4 


-0.040 30103.98(4) 


2 


-0.050 







table lllll (Continued.) 




V = 


V = 1 


V = ^ 


) 




N EnvN 


riy 


/S.E EnvN l^u 


AE EnvN 




AE 


1 16695.884(8) 


8 


-0.014 18308.374(7) 12 


0.006 19857.301(6) 


15 


-0.011 


2 16755.611(11) 


8 


-0.011 18365.854(9) 14 


-0.004 19912.998(10) 


12 


0.002 


3 16844.912(9) 


8 


-0.002 18451.709(7) 13 


-0.019 19996.263(6) 


16 


-0.013 


4 16963.424(13) 


5 


-0.004 18565.166(10) 13 


0.004 20106.670(9) 


14 


0.000 


5 17110.685(9) 


4 


-0.005 18703.797(7) 12 


-0.017 20243.786(7) 


12 


-0.006 


6 17286.063(13) 


4 


-0.013 18884.047(15) 3 


-0.007 20407.009(11) 


7 


0.001 


7 17488.918(10) 


4 


0.002 19076.671(10) 4 


0.009 20595.676(8) 


6 


0.004 


8 17718 453fl4) 


4 


-0.023 19297.282(14) 4 


-0.002 20809.051(12) 


5 


-0.011 


9 17973.850(12) 


3 


-0.030 19543.237(11) 3 


0.023 21046.331(10) 


4 


-0.001 


10 18254.078(18) 


3 


-0.048 


21306.503(13) 


3 


-0.023 


11 18557.648(15) 


2 


-0.078 


21588.485(18) 


1 


0.015 


12 18885.46(2) 


1 


-0.050 


21891.41(3) 


1 


-0.010 


V = 5 


► 


V = 4 








N ErivN 


flu 


/S.E EnvN IT'U 


AE 






1 21343.795(7) 


15 


-0.005 22769.695(7) 14 


-0.015 






2 21397.294(11) 


14 


0.006 22821.398(10) 15 


0.002 






3 21477.206(7) 


16 


0.004 22898.632(7) 15 


-0.012 






4 21583.081(11) 


14 


-0.001 23000.957(10) 16 


-0.007 






5 21714.214(8) 


14 


0.006 23127.895(7) 15 


-0.005 






6 21868.935(14) 


5 


-0.025 23278.629(11) 14 


0.001 






7 22046.649(17) 


1 


-0.009 23451.647(9) 11 


0.003 










jy g 








V = 


V = 1 


V = ^ 


) 




N EnvN 


riy 


AE EnvN fly 


AE EnvN 




AE 


16716.76(3) 


2 


0.040 18225.643(18) 1 


0.017 19649.58(3) 


2 


0.030 


1 16731.26(2) 


5 


0.009 18242.564(13) 6 


0.046 19662.30(2) 


4 


0.070 


2 16761.132(14) 


9 


-0.022 18269.116(9) 9 


0.034 19690.525(19) 


6 


0.045 



table [ml (Continued.) 

V = 3 V = 4 



N EnvN 


riu 


AE EnvN 




AE 


20982.138(17) 


1 


0.052 






1 20996.599(14) 


3 


0.091 






2 21027.241(9) 


6 


0.029 22276.495(14) 


3 


0.085 


3 21075.984(15) 


6 


0.156 22327.149(19) 


2 


0.001 


4 21144.762(10) 


6 


0.038 22397.125(17) 


2 


0.125 


5 21234.849(15) 


4 


0.071 22487.30(2) 


1 


0.079 


6 21347.112(15) 


3 


0.108 22598.58(2) 


1 


-0.010 


7 21481.78(5) 


2 


0.070 






8 21639.07(3) 


1 


0.030 






9 21818.86(3) 


1 


-0.030 






10 22020.93(3) 


1 


0.070 







V = 




V = 1 




V = 2 






N EnvN 


AE EnvN 




AE EnvN ny 


AE 


1 16968.914(17) 


4 


-0.004 18497.025(19) 


5 


0.025 19952.756(14) 


6 


0.034 


2 17017.47(2) 


5 


-0.032 18545.33(2) 


9 


0.076 20000.436(15) 


9 


0.124 


3 17091.044(17) 


5 


-0.014 18618.163(16) 


8 


0.017 20072.116(12) 


8 


0.054 


4 17190.22(2) 


5 


0.007 18715.65(2) 


7 


0.074 20167.775(15) 


7 


0.055 


5 17315.208(16) 


4 


-0.068 18838.130(16) 


7 


0.020 20287.31(2) 


4 


0.190 


6 17466.29(2) 


5 


-0.034 18985.45(2) 


5 


0.024 20430.99(3) 


2 


0.110 


7 17643.211(17) 


4 


-0.091 19157.50(2) 


4 


0.010 20598.47(4) 


2 


0.010 


8 17845.71(2) 


4 


-0.000 19353.99(2) 


4 


0.030 20789.05(7) 


1 


0.060 


9 18073.46(2) 


4 


-0.120 19574.38(2) 


3 


0.010 21002.86(6) 


2 


-0.500 


10 18325.64(3) 


3 


0.010 19818.25(5) 


1 


0.040 






11 18601.48(3) 


3 


-0.140 20084.73(5) 


2 


0.110 






12 18900.26(4) 


3 


0.060 










13 19220.93(6) 


2 


-0.210 










14 19562.52(5) 


1 


0.080 
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table [ml (Continued.) 



V = 




V = 1 




V = 2 






N EnvN riy 


AE' -ErMjAT fly 


AE EnvN n-u 


AE 


1 16998.00(2) 


4 


0.024 18539.50(2) 


5 


0.040 






2 17086.474(17) 


5 


-0.014 18617.862(19) 


6 


0.008 20033.06(2) 


5 


0.070 


3 17202.627(19) 


5 


-0.007 18725.790(19) 


7 


0.010 20130.70(3) 


3 


0.080 


4 17343.627(15) 


6 


-0.047 18860.377(15) 


7 


0.033 20254.70(3) 


5 


0.040 


5 17508.18(2) 


4 


-0.026 19019.62(2) 


7 


0.048 20403.22(4) 


3 


-0.030 


6 17695 856fl5) 


5 


-0 086 19202 388('17) 


7 


-0 018 20597 1A(5) 


1 


0.050 


7 17906.63(2) 


4 


-0.011 19407.97(2) 


6 


-0.044 20795.15(5) 


1 


0.140 


8 18140.58(2) 


3 


-0.110 19635.73(2) 


5 


0.030 20984.82(5) 


1 


0.050 


9 18397.46(2) 


3 


0.010 




21258.76(6) 


1 




10 18676.88(5) 


2 


0.000 










11 18978.67(4) 


1 


0.040 










V = 3 




V = 4 










N EnvN riy 


AE EnvN fly 


AE 






1 21344.541(12) 


2 


0.059 22632.03(5) 


1 


0.030 






2 21408.795(16) 


3 


0.025 22704.67(5) 


1 


0.150 






3 21501.76(2) 


2 


0.096 22791.80(3) 


2 


0.060 






4 21621.287(18) 


3 


0.103 22902.73(5) 


1 


0.010 






5 21765.15(2) 


2 


0.090 23037.01(4) 


1 


0.080 






6 21932.05(3) 


1 


0.060 










7 22120.72(5) 


2 


0.060 














{lsaM5)fA- 








V = 




V = 1 




V = 2 






N EnvN ny 


AE EnvN flu 


AE EnvN n-u 


AE 


2 17401.30(2) 


3 


-0.011 18988.127(15) 


7 


0.023 20511.143(13) 


8 


-0.003 


3 17501.830(13) 


5 


-0.030 19082.632(11) 


9 


0.008 20599.883(9) 


7 


0.047 


4 17633.75(2) 


3 


-0.011 19207.051(15) 


8 


0.039 20717.119(14) 


6 


0.091 


5 17795.867(14) 


4 


-0.117 19360.402(16) 


6 


-0.002 20862.015(15) 


5 


0.085 
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table mi (Continued.) 



{lsa3d6)fAg 



V = 3 




V = 4 




V = 5 




N EnvN fly 


AE EnvN riy 


AE EnvN 


riy AE 


2 21971.433(16) 


4 


0.087 23370.57(3) 


3 


0.070 24709.17(5) 


2 0.070 


3 22054.903(15) 


3 


0.087 23448.57(3) 


3 


0.090 




4 22165 210('16') 


4 


0.080 23551.78(4) 


2 


0.060 




5 22301 619fl5) 


2 


0.041 23679.34(4) 


2 


0.030 




6 22463.39(2) 


2 


0.040 23830.38(4) 


2 


0.060 




7 22649.128(18) 


3 


0.032 








8 22858.24(3) 


3 


-0.020 








9 23089.12(2) 


3 


0.060 












{lsaM5)fA+ 






V = 




V = 1 




V = 2 




N EnvN riy 


AE EnvN ny 


AE EnvN 


riy AE 


2 17401.988(16) 


5 


-0.028 18988.476(15) 


6 


0.034 




3 17504.77(2) 


4 


-0.009 19084.110(17) 


6 


0.031 20600.596(17) 


6 0.084 


4 17641.516(16) 


4 


-0.046 19211.121(14) 


6 


0.009 20719.229(13) 


6 0.051 


5 17811.24(2) 


4 


-0.021 19368.829(19) 


5 


0.041 20866.399(17) 


6 0.051 


6 18012.486(19) 


4 


-0.086 19556.203(17) 


4 


-0.003 21041.275(15) 


6 0.085 


7 18243.37(2) 


4 


0.011 19771.96(2) 


3 


0.060 21242.90(2) 


3 0.135 


8 18502.08(2) 


3 


-0.100 20014.53(2) 


2 


-0.020 21469.95(2) 


3 0.070 


9 18786.62(2) 


3 


0.020 20282.40(6) 


2 


-0.140 21721.25(4) 


2 0.060 


10 19095.39(2) 


3 


-0.130 20573.72(3) 


2 


-0.010 21963.93(6) 


1 0.180 


11 19426.70(4) 


2 


0.070 




22251.87(5) 


2 0.060 


12 19778.99(6) 


2 










14 20542.39(7) 

V = 3 


1 










N EnvN 


AE 








3 22054.98(2) 


1 


0.219 








4 22166.237(14) 


2 


0.073 









